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Specific heat experiments on single crystals of the S—l quasi-one-dimensional bond-alternating 
antiferromagnet Ni(C9H24N4)(N02)C104, alias NTENP, have been performed in magnetic fields 
applied both parallel and perpendicular to the spin chains. We have found for the parallel field 
configuration that the magnetic specific heat (Cmag) is proportional to temperature (T) above a 
critical field He, at which the energy gap vanishes, in a temperature region above that of the long- 
range ordered state. The ratio Cmag/T increases as the magnetic field approaches He from above. 
The data are in good quantitative agreement with the prediction of the c=l conformal field theory in 
conjunction with the velocity of the excitations calculated by a numerical diagonalization, providing 
a conclusive evidence for a Tomonaga-Luttinger liquid. 



PACS numbers: 75.40.Cx, 75.10.Jm, 75.50.Ee 

One-dimensional antiferromagnets and their labora- 
tory counterparts, quasi-one-dimensional (quasi-lD) an- 
tiferromagnets, with an energy gap above a singlet 
ground state are of current interest in condensed matter 
physics. When an external magnetic field is applied to 
any of these systems, one of the excited triplet branches 
goes down in energy due to the Zeeman effect, and the 
energy gap vanishes at a critical magnetic field He- This 
is a quantum critical point at which a transition from a 
gapped spin liquid to a Tomonaga-Luttinger liquid (TLL) 
is expected to occur P, I3, if the system possesses a 
Heisenberg or XY symmetry. The TLL is characterized 
by a gapless A;-linear energy dispersion with an incom- 
mensurate fco and a spin correlation having an algebraic 
decay. 

To date, however, all experimental evidence in favor 
of this scenario has been either controversial or circum- 
stantial. Divergence of the NMR spin-lattice relaxation 
rate and an anomalous hump in the magnetic spe- 
cific heat 0, H, of the spin- 1/2 ladder compound 
Cu2(C5Hi2N2)2Cl4, alias CuHpCl, at low temperatures 
have been interpreted to arise from a TLL. A similar in- 
terpretation of the NMR relaxation rate and the mag- 
netic specific heat has been put forward for the spin- 
1 /2 alternating-bond chain compound pentafuloropheiryl 
nitronyl nitroxide (F5PNN). However, Stone et al. ^Q] 
have proposed on the basis of neutron measurements 
and an examination of the crystal structure that CuH- 
pCl is three-dimensional rather than quasi-lD. Further- 
more, neither compound exhibits a T-linear specific heat 
indicative of the linear, gapless dispersion of a TLL. 

The organometallic Ni-chain compound 
Ni2(C5Hi4N2)2N3(PF6), abbreviated as NDMAP, 



is one of the best-established quasi-lD antiferromagnets 
with an XY symmetry. For this gapped S=l linear-chain 
compound, the H-T phase diagram of the field-induced 
long range order (LRO) in the magnetic fields above He 
for 1 1 chain has been studied in detail and has been 
explained by assuming that the critical exponent of 
the TLL governs the LRO |0. However, no stronger 
evidence for a TLL exists. In this compound, the 
principal axis of the NiNg octahedra which determines 
the crystal-field anisotropy is tilted from the chain axis 
by 15°, with the tilt direction alternating from chain to 
chain. Therefore, neither 7f||chain nor other magnetic 
field directions strictly satisfy an XY symmetry for all 
the chains simultaneously. This lack of axial symmetry 
probably explains why the evidence for a TLL in 
NDMAP is only circumstantial at best. 

A recently synthesized Ni compound 
Ni(C4H24N4)N02(C104), or NTENP ^ for short, 
is an S—1 bond-alternating-chain antiferromagnet. 
NTENP has an energy gap between the singlet ground 
state and the first excited triplet and exhibits an 
LRO above Fc [llEi- Unhkc in NDMAP, the principal 
axis of the Ni^+ ions in NTENP is nearly parallel to 
the chain direction. In this Letter, we report a TLL 
behavior observed in the specific heat of NTENP above 

First, we briefly summarize the crystal and magnetic 
properties of NTENP by referring the reader to Ref. ^3 
for details. NTENP belongs to a triclinic system (space 
group -Pi), and the Ni^+ ions are bridged along the a axis 
by nitrito groups having two alternating bond lengths. 
The Ni chains are well separated by ClO^ anions, thus 
having a good ID nature. The inversion centers are sit- 
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uated in the nitrito groups rather than the Ni^+ ions, so 
that no staggered components of the magnetic moments 
are expected to be retained. The critical fields par- 
allel and perpendicular to the chain direction are 9.3 T 
and 12.4 T, respectively 

The model spin Hamiltionian for NTENP is written as 

'H — y~^[JS2i-i ■ + 5JS2i ■ S2i+i — /ieSigH 

i 

+ D{Stf+E{{Stf-{Syf]l (1) 

where J is the larger of the two exchange constants, 5 
the bond-alternating ratio, g the g tensor of Ni^+, ji-Q the 
Bohr magneton, and D and E the axial and orthorhom- 
bic single-ion anisotropy constants, respectively. From 
the analysis of the magnetic susceptibility data using nu- 
merical results with E=0, the following parameter val- 
ues are known: J/fce = 54.2 K, (5=0.45, D/J ^ 0.25, 
and 11 = 2.14 for the component of g in the chain direc- 
tion ^2|. The ground state of NTENP is expected to 
be a singlet-dimer phase rather than a Haldane phase on 
the basis of the bond-alternating ratio. This expectation 
has been confirmed by magnetic susceptibility, magneti- 
zation, and ESR experiments on Zn-doped NTENP p^ . 
The analysis of the excitations observed by inelastic neu- 
tron scattering in the magnetic fields indicates that the 
strength of E is only 5% of that of D Q. For this rea- 
son, NTENP is an excellent candidate for an observation 
of a TLL. 

We prepared single crystals of hydrogenous and deuter- 
ated NTENP according to the method described in 
Ref.ini The specific heat measurements for deuterated 
and hydrogenous NTENP in the magnetic fields applied 
along the chain direction (a axis) were performed at the 
National High Magnetic Field Laboratory (NHMFL) in 
Tallahassee, Florida. The measurements for a hydroge- 
nous sample in the magnetic fields normal to the chain 
direction were performed at KYOKUGEN, Osaka Uni- 
versity. Some of the results for this field configuration 
have been published in a preliminary report jl3| . The 
lattice contribution of the specific heat was determined 
to be 3.31 xT'^ mJ/K mol by fitting the data at temper- 
atures between 0.2 K and 1.1 K in zero field, where the 
singlet gro und state with a large energy gap of about 
1 meV [l^ makes the magnetic component of the spe- 
cific heat negligible. This lattice contribution has been 
subtracted from the raw data. In addition, we have sub- 
tracted the nuclear contribution of the hydrogen atoms 
from the hydrogenous-sample data at all the fields. 

Figures 1(a) and (b) display the temperature depen- 
dences of the magnetic specific heat for deuterated and 
hydrogenous NTENP samples in the magnetic field ap- 
plied parallel and perpendicular to the chains, respec- 
tively. The hydrogenous sample for chain has given 
similar results to those for the deuterated one. These 
results are not shown here, since they are less extensive 



than the deuterated-sample data. For both field configu- 
rations, a sharp peak signals the LRO above the critical 
field, whose value depends on the field direction. The 
peak position moves up in temperature with an increas- 
ing magnetic field in both cases, but the behavior of the 
magnitude of the peak depends strongly on the field di- 
rection. For 7?||chain, the peak height decreases slightly 
with an increasing magnetic field higher than 12 T as 
shown in Fig. 1(a), whereas it increases steeply with in- 
creasing field for if _L chain (Fig. 1(b)). 

The most striking difference between the results for 
the two field directions is the temperature dependences 
in the region above the ordering temperatures. The spe- 
cific heat for chain is linear in temperature in this 
region, which becomes broad with increasing field. The 
specific heat is only weakly dependent on the magnetic 
field in this region. In contrast, a T-linear specific heat, 
which is indicative of a TLL state, is completely absent 
for iJ_Lchain. The specific heat above the temperature of 
the peak for a given magnetic field decreases much fur- 
ther with increasing temperature than for chain and 
then gradually increases after reaching a broad, rounded 
minimum. A nearly field-independent behavior is found 
at temperatures further away from the peak. It is sig- 
nificant that a T-linear specific heat is not observed for 
_ff_Lchain. In this case, the breaking of the axial sym- 
metry produces a crossover to a 2D Ising transition [l5j , 
causing an absence of a TLL. 

To examine the linear temperature dependence in de- 
tail, we plot in Fig. 2 the specific heat for iJ||chain di- 
vided by temperature for the five designated fields. We 
find that Cmag/T is nearly constant over a temperature 
range above the bottom of the high-temperature side of 
the peak. The upper end of the constant Cmag/T cannot 
be uniquely defined, since this is evidently a crossover 
point instead of a well-defined transition. We take an 
inflection point, which is marked by the arrow, to be the 
upper limit, although this choice is by no means the only 
possibility. The upper boundaries of the TLL marked by 
the arrows in Fig. 2 are plotted in the H-T phase diagram 
of Fig. 3. The diagram clearly indicates the stability of 
the TLL in high magnetic fields, in agreement with the 
theoretical expectation [l^ . 

The values of 7, defined as C^s.g/T in the TLL re- 
gion, are shown in Fig. 4 with solid circles. 7 is weakly 
field-dependent, rising slightly as the magnetic field de- 
creases toward He ~9.3 T. According to the conformal 
field theory for the central charge c=l 7 is given 

by RirkB/Shv, where R and v are the gas constant and 
the velocity of the low-lying excitations, respectively. By 
using a numerical diagonalization method, we have cal- 
culated the dispersion curve of a chain of up to 20 spins 
and have extracted w as a function of the ground-state 
magnetization. The results are sufficiently independent 
of the system size N for 7V=16, 18, and 20. From these 
results and the experimental magnetization curve |l2l |. 
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FIG. 1: (a) Magnetic specific heat of a deuterated sample for 
the magnetic fields applied parallel to the chains, (b) Mag- 
netic specific heat of a hydrogenous sample for the magnetic 
fields applied perpendicular to the chains. 





FIG. 2: Magnetic specific heat of the deuterated sample di- 
vided by temperature for the magnetic fields applied parallel 
to the chains. For clarity, the data have been shifted up by 
50 mJ/K^mol at each field increment. 
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FIG. 3: Magnetic field versus temperature phase diagram of 
NTENP for the magnetic fields parallel to the chains. The 
regions marked LRO and TLL are the long-range ordered and 
Tomonaga-Luttinger liquid phases, respectively. The solid 
and open circles indicate the peak positions of the specific heat 
for the deuterated and hydrogenous samples, respectively, and 
the solid squares show the positions of the arrows in Fig. 2. 
The solid line is the best fit of the deuterated-sample data 
to the expression Tc = A\H - He]" with i/c=9.17 T and 
a— 0.264:. The broken line is a guide to the eye. 



we obtain 7 as a function of the magnetic field. At fields 
up to about 16 T, the agreement with the experimental 
data is excellent as shown in Fig. 4. The calculated re- 
sult at about 20 T is lower than the experimental value 
by about 20%, but the agreement is still impressive. The 
increase of 7 with decreasing H toward He is related to 
the divergence of the density of states at the band edge 
of the spinless fermions in a ID system . The overall 
quantitative agreement between the experimental data 
and the calculation provides a conclusive evidence for a 
TLL. 

The positions of the specific-heat peaks due to the LRO 
are shown in Fig. 3 for _ff ||chain. The peak positions for 
the hydrogenous sample are also plotted for the same 
field direction. As observed in NDMAP the order- 
ing temperature is somewhat lower for the deuterated 
sample, indicating that deuteration slightly weakens the 
interchain exchange interaction. 

The LRO found above He is considered to be the 
Bose-Einstein condensation (BEG) of the triplets [2l| . 
According to the BEG theory |22|, the phase bound- 
ary in the H-T phase diagram obeys the power law 
Tc oc \H — Hcl" . The solid line is the best fit of the data 
for the deuterated sample up to 13 T to this expression 
with 77c=9.17±0.01 and a=0.264±0.002. The previously 
reported a of 0.334 j3| was larger than this, because the 
transition occurred at higher temperature probably due 
to a tilting of the magnetic fields from the chain direction. 
The present Tc values are consistent with those found by 
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FIG. 4: Field dependence of 7=Cmag/r- The solid circles are 
experimental values taken from Fig. 2, and the open circles 
are calculated ones, which are given by RnkB/3hv. Here, R 
and V are the gas constant and the velocity of the low-lying 
excitations, respectively. 



a recent NMR experiment |23|. The a value of NTENP 
is smaller than 0.34 for NDMAP [13 . 

Finally, we comment that the field dependence of the 
peak height of the specific heat for chain is unusual 
from a thermodynamic point of view. Ordinarily, the 
peak at an LRO transition becomes larger as the critical 
region becomes wider with the increasing transition tem- 
perature due to an increasing magnetic field. However, 
the peak height observed for if || chain decreases slightly 
with an increasing magnetic field. The unusual field de- 
pendence probably arises from strong correlations in the 
TLL. In the TLL, the correlation function decays only 
algebraically and the short range order is well developed. 
This leads to a substantial entropy drop above the LRO 
temperature. As the increasing magnetic field makes the 
TLL region wider, the entropy drop becomes larger, leav- 
ing less entropy to be expended at the LRO temperature. 

In conclusion, the specific heat of NTENP provides the 
first conclusive evidence for a TLL in a gapped quasi- ID 
antifcrromagnet. We expect neutron scattering experi- 
ments of NTENP to observe a fc- linear dispersion with an 
incommensurate fco in the TLL phase 0] . But these ex- 
periments present technical challenges, because NTENP 
belongs to a triclinic system, for which the a direction 
(the chain direction) and the direction of the reciprocal 
vector a* are different. 
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